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Abstract 

In the classical one-dimensional solution of fluid dynamics equations all unknown functions depend 
only on time t and Cartesian coordinate x. Although fluid spreads in all directions (velocity vector 
has three components) the whole picture of motion is relatively simple: trajectory of one fluid particle 
from plane x = const completely determines motion of the whole plane. Basing on the symmetry 
analysis of differential equations we propose generalization of this solution allowing movements in 
different directions of fluid particles belonging to plane x = const. At that, all functions but an angle 
determining the direction of particle's motion depend on t and x only, whereas the angle depends on all 
coordinates. In this solution the whole picture of motion superposes from identical trajectories placed 
under different angles in 3D space. Orientations of the trajectories are restricted by a finite relation 
possessing functional arbitrariness. The solution describes three-dimensional nonlinear processes and 
singularities in infinitely conducting plasma, gas or incompressible liquid. 

PACS numbers: 47.10.-g, 47.10.A, 47.65.-d 

AMS classification scheme numbers: 76W05, 76M60, 35C05, 35N10 
Introduction 

The widely-used simplification of fluid dynamics equations is an assumption of one- 
dimensionality of the flow. It is proposed that all unknown functions depend only on 
two variables: time t and Cartesian coordinate x. Motion of fluid particles is allowed 
in all directions, however most of interesting processes (waves of compression and rar- 
efaction, strong and weak discontinuities, etc.) take place along one spatial axis Ox. 
Components of the velocity vector, thermodynamical and all another unknown functions 
are constant on the planes x = const and change from one plane to another. This so- 
lution is often referred to as fluid motion with planar waves. Being comparatively easy 
for an analytical analysis, this simplification provides a great deal of information about 
qualitative properties of fluid motions. However, the classical one-dimensional solution 
can not describe three-dimensional processes in fluid which in fact might be significant 
for the correct description of the picture of the flow. 

In present work we generalize the described classical one-dimensional solution with 
planar waves. In our solution velocity vector is decomposed into two components, one of 
which is parallel and another one is orthogonal to Ox axis. Absolute values (lengthes) of 
the components and both thermodynamical functions (density and pressure) are supposed 
to depend only on t and x. This part of solution coincide with the classical one. However, 
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the angle of rotation of velocity vector about Ox axis is supposed to depend on all inde- 
pendent variables (t,x,y,z). Presence of this function gives the desired generalization of 
the classical solution. 

The proposed representation of the solution was advised by the theory of symmetry 
analysis of differential equations [1, 2]. Indeed, from the symmetry analysis point of view, 
the classical one-dimensional solution is an invariant one of rank 2 with respect to the 
admissible group of shifts along Oy and Oz axis. Whereas the generalized solution is a 
partially invariant one [1] with respect to the full group of plain isometries consisting of 
shifts along Oy and Oz axes and rotation about Ox axis. 

Class of generalized solutions is happened to be a contansive one. It is described by a 
closed system of PDEs with two independent variables, which in the special case coincide 
with classical equations for one-dimensional fluid motions. The angle as a function of four 
independent variables is determined on solutions of the invariant system from a finite (not 
differential) relation, which has a functional arbitrariness. The finite relation allows clear 
geometrical interpretation. This gives opportunity to construct a desired type of fluid 
motion by choosing appropriate arbitrary function in the functional relation. 

Plasma flow governed by the solution possesses a remarkable property. Fluid particles 
belonging to the same initial plane x = const at some moment of time circumscribe the 
same trajectories in 3D space and have identical magnetic field lines attached. However, 
each trajectory and magnetic field line has its own orientation, which depends on the 
position of the fluid particle in the initial plane. The orientation is given by the finite 
relation with functional arbitrariness. Thus, with the same shape of trajectories and 
magnetic field lines one can construct infinitely many pictures of fluid motions by varying 
admissibly directions of particles spreading. 

Intensively studied in recent scientific literature solution of ideal compressible or in- 
compressible fluid equations which is called "singular vortex" or "Ovsyannikov vortex" 
[3, 4, 5, 6, 7, 8, 9, 10] can be treated as the analogous generalization of one-dimensional 
motion with spherical waves. In this solution absolute values of the tangential and nor- 
mal to spheres r = const components of velocity vector field depend only on time t and 
distance r to the origin. An angle of rotation of the vector field about the radial direction 
Or is a function on all independent variables. This solution also allows symmetry inter- 
pretation as the partially invariant one with respect to the admissible group of sphere 
isometries, i.e. of rotations in R 3 . 

The generalized one-dimensional solution with planar waves for ideal gas dynamics 
equations was first obtained in [13]. For all we known, it was not analyzed in details for its 
physical content. In present work we observe equations of ideal magnetohydrodynamics. 
Cases of ideal gas dynamics and ideal liquid can be obtained in limits of zero magnetic 
field H = and constant density p = const respectively. 

The paper is organized as follows. We start from the formulation of the representation 
of solution, which is prescribed by symmetry properties of the main model of ideal mag- 
netohydrodynamics. Substitution of the representation of the solution into the system of 
equations brings a highly-overdetermined system of PDEs for the non-invariant function 
- angle of rotation of the vector fields about Ox axis. Investigation of the overdetermined 
system reveals two main cases, when some auxiliary function h is equal or not equal to 
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zero. From the mechanical point of view these two cases correspond to the compressible 
or incompressible (divergence-free) vector field which is obtained as a projection of the 
velocity field into Oyz plane. In both cases the overdetermined system is reduced to some 
compatible invariant subsystem of PDEs with two independent variables and a finite im- 
plicit relation for the non-invariant function. We give geometrical interpretation of the 
finite relation, which allows keeping track of the singularities, which may take place in 
the flow. We prove that particles trajectories and magnetic field lines are planar curves. 
Moreover, these curves are the same for all particles, which start from the same initial 
plane x = const. This gives opportunity to construct a pattern of the trajectory and 
magnetic field line. The complete 3D picture of the flow is obtained by attaching the 
pattern to every point in fixed Oyz plane in accordance to the directional field defined by 
the finite relation for the non-invariant function. Remarkable, that the same pattern of 
magnetic line and trajectory attached to different directional field in Oyz plane produces 
variety of pictures of plasma motion in 3D space. As an example, the solution is used 
for explicit description of the plasma flow in axisymmetric canal with curved conducting 
walls. 

1 Representation of the solution and preliminary analysis 
1.1 Representation of the solution 

The system of ideal magnetohydrodynamics (tension comes to pressure, thermal conduc- 
tivity is zero, electric conductivity is infinite) has the form [17] 



Here u = (u, v, w) is the fluid velocity vector, H = (H, K, L) is the magnetic vector field; 
p and p are pressure and density. The state equation p = F(S,p) with the entropy S 
gives rise to function A(p, p) defined by A = p(dF/dp). All unknown functions depend 
on time t and Cartesian coordinates x = (x, y, z). 

In the case of arbitrary state equation p = F(S,p) equations (1. 1)— (1.5) admit 11- 
dimensional Lie group Gn of point transformations, which is 10-dimensional Galilean 
group extended by the homothety [11, 12]. Optimal system of subgroups QGn was 
constructed in [15, 16], see also [14]. Examination of QGn shows, that the partially 
invariant solution of described type is generated by 3-dimensional subgroup G3.13 C Gu 
with Lie algebra L 313 spanned by the infinitesimal generators {d y , d z , zd y — yd z + wd v — 
vd w + L&k — K8l} (we use the subgroups numeration according to [14]). 

Indeed, Lie group G3.13 is spanned by shifts along Oy and Oz axes and simultaneous 
rotations about the first coordinate axis in M 3 (x), M 3 (u), and M 3 (H). Invariants of this 
group of transformations in the space of independent variables and dependent functions 



D p + pdiv u = 0, 

D u + p^Vp + p _1 H x rot H = 0, 

Dp + A(p, p) div u = 0, 




divH = 0, D = d t + u-V. 
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Figure 1: Representation of velocity vector u and magnetic field vector H in the partially invariant 
solution. All functions but u> depend on t and x, whereas lo = u>(t, x, y, z). 

M 4 (t,x) x M 8 (u,H,p,p) are 

t, x, u, V = Vv 2 + w 2 , p, p, H, N = VK 2 + L 2 , and vK + wL. (1.6) 

The last invariant may be treated as angle a between the projections of vectors u and H 
into Oyz plane (see figure 1). The general theory of partially invariant solutions may be 
found in [1] . The representation of partially invariant solution is obtained by assigning a 
functional dependence between the group invariants (1.6). In particular, for the solution 
of rank 2 (two invariant independent variables) and defect 1 (one non-invariant function) 
it gives the following representation of solution: 

u=U(t,x), H = H(t,x), 

v — V(t,x) cosuj(t,x,y, z), K — N(t, x) cos (u(t, x, y, z) + a(t, x)j, 

w = V(t, x) sinu;(t, x,y,z), L = N(t, x) sin (u(t, x, y, z) + cr(t, x)) , 

p = p(t,x), p = p(t,x), S = S(t,x). 

Here only the non-invariant function u(t,x,y, z) depends on all original independent 
variables. Functions U, V, H, N, a, p, p are invariant with respect to G3.13- They depend 
only on invariant variables t and x. These functions will be referred to as invariant ones. 
The system of equations for determination of invariant and non-invariant functions will 
be called the submodel of the main model of ideal magnetohydrodynamics. 

1.2 Analysis of the submodel 

Substitution of the representation (1.7) into (1. 1)— (1.5) gives the following result. The 
continuity equation (1.1) allows introduction of new unknown invariant function h(t,x), 
defined by the following relation 

Dp + p(U x + hV) = 0. (1.8) 
Hereinafter D denotes the invariant part of the differentiation along the trajectory 

D = d t + ud x . 
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The remaining part of the continuity equation gives an equation for function uj: 

sincu u)y — coscj uj z + h = 0. (1.9) 

Another equations for invariant functions follow from the first components of momentum 
(1.2) and induction equations (1.4), and also pressure equation (1.3). 

DU + p- 1 Px + p- 1 NN x = 0, (1.10) 
DH + hHV = 0, (1.11) 
Dp + A(p,p)(U x + hV) = 0. (1.12) 

The rest of system (1.1)— (1.5) gives rise to the overdetermined system for function uj. 
From a nondegenerate linear combination of equations (1.2) in projections to Oy and Oz 
axes one obtains 

pVuj t + (p UV — HN cos a) uo x + {pV 2 cos uj — N 2 cos a cos(uj + cr)) uo y (1-13) 
+ (pV 2 sin a; — iV 2 cos cr sin(w + cr)) uj z — H(N X sin cr + N cos crcr x ) = 0. 

HN sin cr + iV 2 sin cr cos(cj + cr) + iV 2 sin a sin(cj + a) cj z (1-14) 
+p 5\/ + i^iV sin aa x - HN X cos cr = 0. 
The same operation with remaining two induction equations (1.4) provides 

Nu t + (NU - HVcosa)uj x + VN sin cr sin(cu + cr) ^ (1.15) 
-VJV sin cr cos(cj + cr) w 2 + OTcr + HV X sin cr = 0. 

sin cr c<j x + A 7 "!/ cos cr sin (c<j + cr) uo y (1-16) 
-A^V cos cr cos(cj + cr) cj 2 - DN + i7\4 cos cr - NU X = 0. 
Finally, equation (1.5) is transformed to 

A(sin(u; + cr)uJy — cos(c<j + cr) lu z ) — H x = 0. (1-17) 

The overdetermined system (1.9), (1.13)— (1.17) for non-invariant function uj should 
be investigated for compatibility [18]. At that we observe only solution with functional 
arbitrariness in determination of function uj. This condition, in particular, implies non- 
reducibility of the solution to the classical one-dimensional solution with planar waves. 

Function uj determines with only constant arbitrariness if it is possible to express all 
first-order derivatives of uj from the system of equations (1.9), (1.13)— (1.17). In order to 
prohibit this situation one should calculate a matrix of coefficients of the derivatives of 
function uj and vanish all its rank minors. This leads to the following four cases: 

l.H = 0; 2. A = 0; 3. V = 0; 4. a = or cr = tt. (1.18) 
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By definition (1.7) functions V and iV are non-negative. Values a = n and a = in the 
case 4 (1.18) differ only by the sign of function N. Both can be observed in the same 
framework for a = 0, non- negative V and arbitrary N. 

Cases 2 and 3 in classification (1.18) correspond to the magnetic field or velocity 
parallel to Ox- axis. Both of them embed into the case a = 0. Indeed, if a = then the 
velocity vector at each particle and its magnetic field vector belong to the plane, which 
is orthogonal to Oyz coordinate plane. Therefore, cases 2 and 3 are degenerate versions 
of this more general situation. Case 4 will be observed as the main case in the following 
calculations. In case of pure gas dynamics H = three of four conditions (1.18) satisfied 
automatically, hence the solution is irreducible without any additional restrictions. 

1.3 Case of planar magnetic field 

Let us first observe the case H = 0, when the magnetic field vector is parallel to Oyz 
plane. The compatibility condition of equations (1.9) and (1.17) in this case is 

( cos(cj + er) uj y + sin(w + a) uj z ) h = 0. (1-19) 

For h — the determinant of the homogenous system of algebraic equations (1.9), (1-17) 
for uj y and uj z is sin a. Hence, the solution is non-trivial only for er = or a = ir. The 
case h ^ leads to the reduction following from equations (1.17) and (1.19). Thus, the 
non-trivial solution exists only for sin a = 0, i.e. case 1 in the classification (1.18) contains 
in case 4. 



2 The main case /i / 
2.1 Equations of the submodel 

Let us observe the main case a = 0. From the mechanical point of view it corresponds 
to a plasma flow where velocity and magnetic field vectors at each particle are coplanar 
to Ox axis. Equations (1.8), (1.10)-(3.36) belong to the invariant part of the submodel. 
Besides, equation (1.14) gives 

D V - p^HNx = 0. (2.20) 
From equation (1.16) taking into account (1.9) one obtains 

DN + NU X - HV X + hNV = 0. (2.21) 

Finally, equation (1.17) due to the relation (1.9) can be written as 

H x + hN = 0. (2.22) 

In addition to the equation (1.9), the non-invariant part of the determining system con- 
tains two equations, which follow from (1.13), (1.15): 

pVu t + (pUV - HN) uj x + (pV 2 - N 2 ) (cos ww y + smuw z ) = 0, (2.23) 
Nu t + (NU - HV) u x = 0. (2.24) 
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Elimination of the derivative oj t from equations (2.23), (2.24) gives the following classifying 
relation 

{pV 2 - N 2 ) (Hu x + iV(cos u uj y + sin u u z )) = 0. (2.25) 

We observe only the case when the second factor in (2.25) vanishes. The compatibility 
conditions of equations (1.9), (2.24), and (2.25) are 

ND h — HVh x = 0, (2.26) 
Hh x + h 2 N = 0. (2.27) 

For h 7^ there is an integral 

H = H h, H = const. (2.28) 
Thus, the submodel's equations are reduced to the following ones. 

Dp + p(U x + hV) = 0. (2.29) 

D U + p~ l p x + p~ l NN x = 0. (2.30) 

DV- p- 1 H o hN x = 0, (2.31) 

Dp + A(p,p)(U x + hV) = 0, (2.32) 

DN + NU X - H hV x + hNV = 0, (2.33) 

D h + Vh 2 = 0, H h x + hN = 0. (2.34) 

The obtained system (2.29)-(2.34) inherits the overdetermination of the initial MHD 
equations (1.1)— (1.5). However, its compatibility conditions satisfied by virtue of the 
system itself. Indeed, the only nontrivial compatibility condition of the system (2.29)- 
(2.34) is given by two equations (2.34) for function h. Cross-differentiation of (2.34) 
shows that their compatibility condition coincide with equation (2.33), i.e. is already 
contained in the system. The most general Cauchy problem for system (2.29)-(2.34) 
requires assigning functions p, U, V, p, N at t = as functions of x, and fixing a constant 
value of h at t — over some plane x = const. For pure gas dynamics H = the second 
equation (2.34) satisfies identically, hence the initial data for h become h(0,x) = h (x). 
System (2.29)-(2.34) equipped by the suitable initial data can be solved numerically. It 
also allows exact reductions to systems of ODEs since the admitted symmetry group is 
obviously nontrivial. 

Equations (2.23)-(2.25) for the non-invariant function can be integrated. Function uo 
determines by the following implicit equation 

F(y — t cos ou, z — t sin oS) = (2.35) 

with t = 1/h and arbitrary smooth function F. In case of pure gas dynamics H = 
equation (2.24) identically satisfied. Therefore, function F in the general formula (2.35) 
for function uj also arbitrarily depends on ^: £ t + U£ x — 0. Results of the performed 
calculations are summarized in the following theorem. 

Theorem 2. In the main case a = and h ^ the invariant functions are determined 
by the system of differential equations (2.29)-(2.34)- The non-invariant function uj is 
given by the implicit equation (2.35) with arbitrary smooth function F. 
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Figure 2: Geometric interpretation of the solution ui — u(r(t,x),y, z) of the implicit equation (2.35). 
Curve 7 : F(y, z) = is determined by the same function F as in (2.35). Function u> at given point M 
is the angle between the direction of line segment AM and Oy axis, where A € 7 and \AM\ = r. 




Figure 3: The field of directions is defined by the implicit equation (2.35) in the stripe of determinacy of 
width 2r with curve 7 : F(y,z) = as a medial line. In this example F = z — smy. At the points of 
limiting equidistants the field of directions is orthogonal to the equidistants. 



2.2 Geometrical construction of the field of directions 

Here we give an algorithm for solving the implicit relation (2.35) over some fixed plane 
x = xo at time t = to- Suppose that function F in (2.35) is fixed. This specifies a 
curve 7 = {(y,z) | F(y,z) = 0}. In order to find angle uj at arbitrary point M = (y,z) 
one should draw a line segment AM of the length r such that A 6 7. The direction 
of AM gives the required angle u as it is shown in figure 2. Function u is only defined 
at points located within distance r from the curve 7. The rest of Oyz plane does not 
belong to the domain of uj. Boundaries of the domain of u are r-equidistants to 7. As 
x grows, function r changes according to the solution of equations (2.29)-(2.34). This 
prescribes modification of the cu-domain over different planes x = const. Thus, the domain 
of function uj (hence, of the whole solution (1.7)) over each plane x = const is a stripe 
of determinacy of the width 2r with curve 7 as a centerline (see figure 3). The stripe of 
determinacy is bounded by equidistants curves to 7. Over the boundaries of cj-domain 
the field of directions uj is orthogonal to the boundaries. 

Inside its domain function uj is multiply-defined. Indeed, there are could be several 
line segments AM with A 6 7 giving rise to several branches of function uj. However, it 
is always possible to choose a single-valued and continuous branch of uj. 
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Discontinuities of iv may appear in cases when the equidistants to 7 have the dovetail 
singularities. The observations illustrated by figure 4 show that every branch of function 
iv necessary have a line of discontinuity inside or at the border of the dovetail. In figure 
4 the curve 7 is a sinusoid shown at the bottom of figures; the curve on the top is the 
equidistant shifted at large enough distance r. For the convenience we draw the circle 
of radius r with center at chosen point M. Each intersection of the circle with 7 gives 
rise to a branch of iv. Let us take M outside of the dovetail (figure a). There are two 
branches of iv at M. As M moves towards the borders of the dovetail, both branches 
change continuously (figure b). At the border of the dovetail the new branch of iv appears 
(figure c). The latter splits into two branches inside the dovetail (figure d). As M reaches 
the right boundary of the dovetail the two "old" branches of iv sticks together (figure e) 
and disappear as M leaves the dovetail (figure f). One can chase, that the branches of iv 
obtained on the right-hand side of the dovetail are different from the ones existed on the 
left-hand side of the dovetail. 

The dovetails do not appear if r < mini?(x), where -R(x) is a curvature radius of 

curve 7 at x. So, one can avoid the singularities either by choosing the solution with 
small enough r or by fixing the curve 7 with large curvature radius. Described disconti- 
nuities takes the solution out of class (1.7). They can not be interpreted in shock waves 
framework. Indeed, over the line of discontinuity only the direction of the magnetic and 
velocity vector fields change, while their absolute values together with thermodynamics 
functions remain continuous. Another type of transverse or alfveic waves [17, 19] char- 
acteristic to ideal MHD equations also can not explain the discontinuity as long as the 
magnetic and velocity fields rotates not across the front of discontinuity. 

Appearance of the dovetail singularities physically mean magnetic field lines, which 
pass through different point in some initial plane x = const collide in their further devel- 
opment. This happens if the function r increases along the magnetic lines such that the 
r-equidistants to 7 became non-smooth. In the vicinity of the collision point the solution 
leaves the prescribed class (1.7); the corresponding fluid flow should be observed either in 
general 3D framework, or in terms of an extended main model, i.e. taking into account 
magnetic or kinematic viscosity as it is observed in magnetic reconnection problems [20]. 
This nonlinear process is specific to the constructed solution, and can not take place in the 
classical one-dimensional solution with planar waves, where all magnetic lines are parallel 
to each other. 

3 Case h = 

3.1 Equations of the submodel 

From the mechanical point of view this case means that the projection of vector field u 
into the plane x = const is incompressible, i.e. its divergence is zero. This case is observed 
separately because the non-invariant function tv is determined by different algorithm. 
For h — integral (2.28) is not valid. Instead, equations (1.11) and (2.22) give 

H = H = const. 
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Figure 4: The behaviour of function uj over the dovetail. There are two branches of u> outside the dovetail 
in Figures (a), (b) and (f); three branches of u> at the borders of the dovetail in Figures (c) and (e); and 
four branches of solution inside the dovetail in Figure (d). 

Thus, equations of the invariant system are 

Dp + pU x = 0, 

DU + p- 1 Px + p~ 1 NN x = 0, 

DV — p~ x HqN x — 0, (3.36) 
Dp + A{p,p)U x = 0, 
DN + NU x -H V x = 0. 

This system of 5 equations serves for determination of 5 unknown functions U, V, N, p, 
and p. The non-invariant function uj is restricted by equations (1.9), (2.24), and (2.25). 
Suppose that its solution uj = uj(t, x, y, z) for N ^ and pV 2 —N 2 ^ is given implicitly by 
the equation $(t, x, y, z, uj) = 0, $ u ^ 0. The system (1.9), (2.24), and (2.25) transforms 
as follows 

$ fe = 0, $ t + U$ x + V$ j = 0, H o ^ x + N^j = 0. (3.37) 
Here Ojk is a Cartesian frame of reference rotated on angle uj about the origin. 

j = ycosuj + zsinuj, k = — ysmuj + z cos uj. (3.38) 

Integrals of system (3.37) are uj and j — <p(t, x), where function <p(t, x) satisfies the overde- 
termined system 

(fH + U(p x = V, H <p x = N. (3.39) 

The compatibility condition of equations (3.39) is the last equation of the invariant system 
(3.36). Differential one-form 

H dip = (H V - NU)dt + Ndx 
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Figure 5: Given a value of ui at some point M, the 
auxiliary Ojk frame of reference is defined as shown. 
The projection of M into the Oj axis is called the 
base point for M. The set of all the base points 
for different M with different lo(M) forms the basic 
curve 7. 



Figure 6: Given curve 7 one can find uj at any point 
M of the corresponding Oyz plane. Circle Sm with 
diameter OM should be drawn. Let Ai be points of 
intersection of Sm with 7. For each Ai the angle uj 
at M is given by the direction OAi as shown. 



is closed, therefore function <p can be found by integration as 

(t,x) 

ip(t, x) = J dip. 

Note, that the initial data for function ip is given by only one constant ip(t Q ,xo). The 
non-invariant function uj can be taken in the form of the finite implicit equation 

j = f(u) + <p(t,x) (3.40) 

with arbitrary smooth function /. The result is formulated in the following theorem. 

Theorem 3. In the case a = h = the invariant functions are determined from 
equations (3.36), (3.39). Function uj is given by the implicit equation (3.4-0). 

3.2 Construction and properties of the field of directions 

Now we clarify a geometrical interpretation of the implicit relation (3.40). Let us fix a 
plane x = xq and time t = to. For simplicity we assume <p(to,xo) = 0. Let the value 
of uj satisfying (3.40) is known at some point M = (y, z) of the plane x = xq. Consider 
a Cartesian frame of reference Ojk turned counterclockwise on angle uj with respect to 
Oyz (see figure 5). By the construction, j-coordinate of point M and angle uj are related 
by j = f{uj). All points with the same coordinate j and arbitrary coordinate k satisfy 
the same relation. A point satisfying the relation (3.40) with zero coordinate k will be 
referred to as the base point for chosen values of j and uj. The locus of all base points for 
various j and uj gives the basic curve 7. On the plane Oyz the basic curve 7 is defined in 
polar coordinates y = rcosO, z = rsin# by the equation r = f(0). Note, that since the 
value of j can have arbitrary sign, both positive and negative values of polar coordinate 
r are allowed in the construction of 7. 
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Figure 7: Curve a 7 is denned by equation r = cos 29 
where both positive and negative values of r are al- 
lowed. Point B corresponds to the part of the curve 
with negative r. The direction V2 assigned to B is 
therefore opposite to the one given by the segment 
OB. 



Figure 8: Point M belongs to the boundary of the 
domain of function lu if the circle Sm is tangent to 
7 at some point A. From the elementary geometry 
vectors OA and AM are orthogonal. This allows 
expressing vector m in terms of x and x x . 



The obtained geometrical interpretation provides an algorithm of construction of the 
vector field, which is defined by the angle uj of deviation from the positive direction of the 
Oy axis. Angle uj is determined from the solutions of implicit equation (3.40). Suppose, 
that function / in equation (3.40) is given. This means, that one can construct the basic 
curve 7 by the formula r = f(8) in polar frame of reference on Oyz plane. Determination 
of angle uj at the point M = (y, z) of the plane x = xq requires the following operations 
as illustrated in figure 6. 

1. Draw a circle Sm with diameter OM. 

2. Find the intersection points A4 of the circle Sm with curve 7. If Sm does not intersect 
7 then M does not belong to the domain of uj. 

3. The angle between the line segment OA t and a positive direction of Ox axis gives a 
value of the angle uj at point M (see figure 6). 

4. Angle uj has the same value at all points of the line passing through the line segment 
AiM. 

As mentioned before, function / can be both positive and negative. Negative / cor- 
responds to negative coordinate j. Thus, if the point of intersection of the curve 7 and 
auxiliary circle Sm belongs to those part of the curve, which corresponds to the negative 
values of /, then the vector field should be taken with the negative sign, i.e. instead 
of uj one should take uj + tt. This situation is shown in figure 7. Here the curve 7 is 
determined by the equation r = cos 26. The circle with diameter OM for M = (3, 3) has 
two points of intersection with curve 7. Point A belongs to the "positive" part of curve 
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Figure 9: The vector field denned by the curve 7 : y 2 + z 2 = R 2 . 

7, therefore it defines the direction vi, codirectional with the segment OA. Point B lies 
on the "negative" part of 7, i.e. the corresponding direction v 2 , is opposite to the one, 
defined by the segment OB. 

Next, it is necessary to find the domain of function u = u(t ,x ,y, z) defined by the 
implicit equation (3.40) over the plane x = Xq. Assume that curve 7 is given. Point M 
belongs to the boundary of the domain if the circle Sm with diameter OM is tangent 
to curve 7 at some point A (see figure 8). Let the position vector of point M be m. 
Parametrization of 7 is taken in the form x = x(s) with some parameter s G A C R. 
From the elementary geometry m = x + ax 1 , where x -1 • x = 0. The tangency condition 
of the circle and curve 7 gives (m/2 — x) ■ x = 0. Here and further the upper dot 
denotes the differentiation with respect to s. Substitution of the expression for m form 
the first equality into the second one provides (ax J -/2 — x/2) ■ x = 0. The scalar a is then 
determined by 

x • x 

a = — j r- 

x^ • x 

Thus, the border of the domain of function uj has the following parametrization 

m = x+ * L ' X x x . x = x(s), sGAcl. (3.41) 

Note, that m does not depend on the choice of the sign and length of x -1 . At the border's 
points the vector field defined by u has x direction. This direction is orthogonal to the 
border. Indeed, 

(I I \ • " • I 

x + dx + ax ) ■ x = x ■ x H : x ■ x = 0. 
x- 1 - ■ x 

The last expression vanishes because from x ■ x -1 = it follows x • x -1 = — x ■ x -1 . 

As an example, let us take 7 to be the circle y 2 + z 2 = R 2 . The border of the domain 
of u in this case coincide with the circle 7 because for each point x of the border one has 
x ■ x = 0. The corresponding vector field describes a flow from the cylindrical source and 
shown in figure 9. In limit R = one obtains a vector field corresponding to the rotation 
around the origin. 
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Figure 10: Trajectories and magnetic field lines are planar curves, which are the same for all particles, 
belonging to the same plane x = const. In order to determine the flow in the whole space it is required 
to set up an admissible vector field of directions in some plane x = xo (i.e. to determine function ui 
consistent with equations (2.35) or (3.40)) and calculate trajectory and magnetic field line for arbitrary 
particle in this plane. The whole picture of the flow is obtained by attaching the trajectory and the 
magnetic line pattern to each point on the plane x = xq in accordance with the vector field of directions. 

4 Particles trajectories and magnetic field lines 
4.1 Trajectory and magnetic field line pattern 

First of all, let us notice that from equations (2.23), (2.24) for pV 2 — N 2 ^ follows the 
equality 

Du = 0. (4.42) 

The trajectory of each particle is a planar curve. Indeed, equation (4.42) implies that 
angle u has constant value along each trajectory. Hence, the whole trajectory belongs 
to the plane, which is parallel to Ox axis and turned on angle uo about this axis. The 
same holds for a magnetic field line, because vanishing of the second factor in (2.25) 
is equivalent to constancy of u along each magnetic curve. Thus, for each particle its 
trajectory and magnetic field line are planar curves, which lie in the same plane defined 
by the angle u>. 

The second important property follows from the representation of the solution (1.7). 
Let us set up a Cauchy problems for trajectory of some particle. The particle moves in 
its plane, hence in this plane the motion is completely defined by components of velocity 
U and V. These two functions depend only on invariant variables t and x. Therefore, 
for any two particles, which belong to the same plane x = xq at initial time t — to the 
Cauchy problems for the trajectories coincide. Of course, the two different particles move 
in their own planes, but both trajectories as planar curves are exactly the same. Similar 
observation is true for any two magnetic lines passing through two different points in the 
same plane x = x . Thus, one can construct a pattern by calculating the trajectory and 
the magnetic field line for any particle in the plane x = xq. The pattern attached to 
each points in the plane x = Xq inside of the domain of function u according to the field 
of directions defined by function uj produces the 3D picture of trajectories and magnetic 
field lines in the whole space. The described algorithm is illustrated in figure 10. 
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In order to construct the pattern let us observe a plane of motion of some particle, 
which is located at initial time t — to at some point M = (xq, yo, z$). This plane is 
parallel to Ox axis and turned about Ox axis on angle u. Cartesian frame of reference 
is defined in the plane of motion as follows. The origin O' of the frame is placed at the 
projection of point M into Oyz plane. One of the coordinate axes is chosen to be parallel 
to Ox axis and denotes by the same letter x. Another axis O'l is placed orthogonally to 
O'x such that the frame O'xl has right orientation (see figure 10). Particle's trajectory 
in this frame of reference is defined by the solution of the Cauchy problem 

dx 

— = U{t,x), x(t )=x . (4.43) 

The dependence x = x(t,xo) given by a solution of (4.43) allows finding the dependence 
I — l(t) along the trajectory by the formula 

t 

l(t) = J V(t,x(t,x ))dt. (4.44) 
to 

The planar curve determined by the dependencies x — x(t, xo) and I — l(t) forms a pattern 
of the trajectory for any particle, which belongs to the plane x = x at t = t . Equations 
of particle's trajectory in initial Oxyz-fr&me are restored in the form 

x — x(t, xo), y = yo + lit) cosc^o, z = z + l(t) smuo- (4.45) 

Here uo = uj(to, xo) is the value of angle u taken at initial time t = to at point M. 

The magnetic field line at t — to is an integral curve of the magnetic vector field. The 
pattern of the magnetic curve passing at t — t through the plane x = x is given by 

l(x)= [^Ads 

xo 

Equations of the magnetic field curve in Oxyz frame of reference are restored as 

f N (to,s) , . f N(t ,s) , 

y = y + cosuo / — — -ds, z = z + sin w / 777: r«S- (4.46) 

J H(t ,s) J H(t ,s) 

xo xo 

Derivation of these formulae is similar to those given for the trajectory (4.45). 
Thus, the following properties of plasma motion holds (see figure 10). 

• Trajectories and magnetic lines lie entirely in planes, which are orthogonal to the 
Oyz-p\ane and turned on angle u about Ox axis. 

• All particles, which belong at some moment of time t — to to a plane x = xq, 
circumscribe the same trajectories in planes of each particle motion. Magnetic field 
lines passing through a plane x = Xo are also the same planar curves. 

• Angle of rotation about Ox-axis of the plane containing the trajectory and the mag- 
netic line of each particle is given by function u>, which satisfies equation (2.35) or 
(3.40). 
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4.2 Domain of the solution in 3D space 



The constructions above show that the whole area in 3D space occupied by the solution 
is obtained as follows. In fixed plane x = x function u has some (in many cases, finite) 
definition domain, bounded by r-equidistants to 7 for h 7^ and by the curve (3.41) for 
h = 0. In both cases the field of direction defined by u in x = Xq plane is orthogonal to 
the boundary of the cu-domain. In order to obtain boundaries of the whole 3D domain of 
the solution one should attach the magnetic line pattern, calculated on some particular 
solution of the invariant system, to every point of the boundary of cu-domain in plane 
x = x according to the usual algorithm. This gives a canal woven from magnetic lines 
which pass through boundaries of the u(tx , y, z)-domain and intersect x = x plane. 
The walls of the canal can be interpreted as rigid infinitely conducting pistons. Due to 
the well-known property of magnetic field line freezing-in, the walls are impermeable for 
plasma. In case of stationary solution the walls are fixed. In non-stationary case the walls 
extend or shrink according to the behavior of function r for h 7^ and ip for h = 0. In 
case of finite cu-domain (it can always be restricted to a finite one) each rr-cross-section of 
the 3D-domain of the solution is finite, therefore both magnetic and kinetic energy have 
finite value in each rc-layer. 



4.3 Stationary flow 

As an example we observe a stationary solution of system (2.29)-(2.34). Suppose that all 
unknown functions depend on x only. This leads to the following system of ODEs: 

Up' + p(U' + hV) = 0. (4.47) 

UU' + p- l p' + p~ l NN' = 0. (4.48) 

UV' - p^HohN' = 0, (4.49) 

Up + A(p, p) (U' + hV) = 0, (4.50) 

UN' + NU' - H hV + hNV = 0, (4.51) 

Uti + Vh 2 = 0, H h' + hN = 0. (4.52) 

Elimination of the derivative h' in equations (4.52) gives the finite relation 

H Vh = UN, (4.53) 

which states collinearity of the magnetic and velocity fields at each particle. The same 
property holds for the analogous spherical solution [9]. Equation (4.51) is satisfied iden- 
tically by virtue of (4.53). 

Equation (4.50) gives entropy conservation 

S = S . (4.54) 

Equation (4.47) under condition (4.53) gives the flow rate integral 

pU = nh, n = const. (4.55) 

Substitution of the obtained integrals into (4.49) allows finding the following relation 
between the tangential components of velocity and magnetic fields 

nV — H N = m, m = const. (4.56) 
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Integration of equation (4.48) gives the Bernoulli integral 

U 2 + V 2 + 2 / ^ = b 2 , b = const. (4.57) 
J P 

The only equation left to integrate is any of two equations (4.52). With its aid all 
unknown functions may be expressed in terms of the "potential" r = 1/h as 

^ mT + g ° V , V = "^ + g ° V , iV=^, " V , . (4.58) 

htt' nr r r mr + if^r' 

a) Let m ^ 0. Using the admissible dilatations it is convenient to make m = n = 
sign(rr'). Expressions (4.58) become 



T + H 2 T> T + H 2 T> H H T> T> 

U=—^, V=—, H=- N=— P = 7TW - (4.59) 



Substitution of (4.59) into the Bernoulli integral (4.57) produces an equation for r. In 
case of polytropic gas with the state equation p = Sp 1 it has the following form 

This ODE for t(x) is not resolved with respect to the derivative r', which complicates its 
investigation. Examples of analysis of such non-resolved ODEs can be found in [4]-[7]. 
One can show that there are several branches of solution t(x) of equation (4.60) passing 
through each point in (x, r) plane, which correspond to different relations between the 
velocity U and the characteristics speeds of MHD system (1.1)— (1.5). It is possible to 
switch between different branches of the solution via fast or slow shock waves. However, 
this investigation lies outside of the scope of this paper. 

b) In case m = after some straightforward simplifications we obtain the following 
solution of system (2.29)-(2.34): 

U = Hq sech x, V = Hgtanhx, r = cosh a;, 

(4.61) 

H = H sech x, N = H tanhx, p = H ~ 2 , S = S . 

One can check that (4.61) represents a special case of the more general S. Chandrasekhar 
solution [21]. This solution is also invariant with respect to infinite group of Bogoyavlenskij 
transformations [22]. The simplicity of solution (4.61) gives opportunity to use it for 
demonstration of geometrical algorithms given in previous sections. 

Streamlines and magnetic field lines coincide and are given by formulas (4.45) with 
xq = and 

l(x) = coshx - 1. (4.62) 

In each plane of particle's motion the streamline is a half of catenary. Note, that solution 
(4.61) can be continuously adjoined with the uniform flow along Ox axis. Indeed, in 
section x — all functions in (4.61) and their derivatives take values compatible with 
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Figure 11: Field of direction obtained by the algorithm of section 2.2. Here 7 is the middle circle of 
radius R. Three cases according to the relation between R and r distinguishes. In all cases the domain 
of the vector fields is an annular stripe of determinacy between two equidistant curves (inner and outer 
circles in the diagrams). 



the uniform flow. Let us construct a solution, which switches the uniform flow to the 
generalized one-dimensional solution (4.61) at the section x = 0. The corresponding 
streamline is a straight lines for x < and a half of the catenary for x > 0. In order 
to get the whole three-dimensional picture of motion this streamline pattern should be 
attached to each point of the plane x = according to the direction field defined by 
function uj. 

Function uj is determined by the implicit equation (2.35). Algorithm of section 2.2 
requires assigning some particular function F, or some curve 7 : F(y, z) = 0. Let the 
curve 7 be a circle y 2 + z 2 = R 2 . The corresponding function u is determined at each 
point of the plane x = by equation (2.35). Figure 11 shows the vector fields obtained for 
different relations between r and R. For R > t the vector field is defined in the annular 
area between two circles of radii R ± r. On the inner equidistant circle |x| = R — r the 
vector field is directed outside of the stripe of determinacy towards the origin. In case 
R = t the inner equidistant circle shrinks into the origin x = 0. At that, the vector field 
becomes multiply-determined at this point. Finally, for R < r the inner equidistant turns 
inside out and becomes a circle of radius r — R with the vector field on it directed inside 
of the stripe of the determinacy. These three vector fields generate different pictures of 
motion in whole 3D space. 

The streamline pattern described above should be attached to each points of Oyz plane 
inside the corresponding domain of u according to the directional fields shown in figure 
11. Because of the obvious central symmetry of the vector fields the whole picture of 
motion is axially-symmetrical. The axial section of the area in 3D space, occupied by the 
corresponding flow is shown in figure 12. 

We assume that uniform flow for x < changes at section x = to the flow, described 
by the solution (4.61). Depending on the relation between r(0) and R three different 
pictures of motion are possible. Each particle moves along the same planar curve, however 
orientation of the streamlines in the space differ from one particle to another. Three- 
dimensional visualization of the motion for R > r(0) is shown in figure 13. 
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Figure 12: Axial sections of axially-symmetrical canal occupied by the plasma flows. The uniform flow 
in cylindrical canal for x < switches at section x — to the flow in the curvilinear canal for x > 
described by the solution (4.61). The boundary of the canal is a rigid wall. Cases a, b and c correspond 
to the vector fields in figure 1 1 . In the diagrams a and c the canal has an inner cylindrical core. 

Conclusion 

In present work a new solution of ideal fluid dynamics equations, describing three-dimen- 
sional motions of plasma, gas and liquid is constructed. The solution is determined by 
a system of equations with two independent variables, which is analogous to the classi- 
cal system for one-dimensional fluid motions. At that, the new solution describes spatial 
nonlinear processes and singularities, which are impossible to obtain in the classical frame- 
work. 

In the constructed solution particles trajectories and magnetic field lines are flat curves. 
Trajectory of each curve and its magnetic field line belong to the same plane parallel to 
Ox axis. In contrast to the classical one- dimensional solution, plane of motion of each 
particle has its own orientation, which is given by an additional finite relation. The 
functional arbitrariness of the finite relation allows varying the geometry of obtained 
motion in accordance to the problem under consideration. Depending on the chosen 
geometry, singularities on the border of the region, occupied by fluid, may appear. In 
such cases particles may collide at the border of the domain of the flow. The criterion of 
singularities appearance in terms of invariant properties of the arbitrary function, which 
specifies the geometry of the flow is given. 

The obtained solution may be used as a test for numerical modeling of complicated 
three-dimensional flows of infinitely conducting plasma. It also may serve for theoretical 
investigations of three-dimensional singularities of the ideal fluid and plasma motions. 
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Figure 13: Tree-dimensional visualization of motion. Fragments of the canal's walls and the streamlines 
are shown. Each streamline has a shape of the same flat curve. Orientation of each streamline is defined 
by the vector field in figure 11a. The axial section of the canal is represented in figure 12a. 
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